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Overview 


Stennis Space Center 



• Objective 


Evaluate accuracy of a prototype algorithm that uses satellite- 
derived atmospheric products to generate scene reflectance 
maps for high spatial resolution (HSR) systems 


• Approach 

- Implement algorithm in an end^to-end process 

- Compare algorithm generated scene reflectance maps with 
ground-truth data 

- Identify algorithm sensitivities 
Provide recommendations 


• Constraints 

Ground truth available only in VNIR spectral range 
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Atmospheric Correction 

Stennis Space Center 

• Atmospheric correction is the process of converting satellite signals 
(at-sensor radiance) to ground reflectances 

- Removes atmospheric and solar illumination effects 

• Benefits 

Improves change detection 

- Used with spectral library based classifiers 

- Simplifies satellite data intercomparisons 

• Different levels of atmospheric correction yield different 
approximations of scene reflectance 

- Planetary reflectance - no knowledge of atmosphere 

- Ground reflectance using knowledge of atmosphere 

- Ground reflectance using knowledge of atmosphere and adjacency 
effects 
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Planetary Reflectance 

— ^ ^ ■— ^ Stennis Space Center 

First-order approximation - no knowledge of atmosphere 



Where : 

p p = Planetary reflectance 

L T0A = Top of atmosphere (at -sensor) radiance 

6 = Solar zenith angle 

E sun = Solar exoatmospheric irradiance 

d = Sun - Earth distance 
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Atmospheric Correction Algorithm 

Implementations 

Stennis Space Center 

• Use knowledge of atmosphere to determine the 

constants necessary to convert satellite signals to scene 
reflectances 

- Ground-based reflectance measurements (direct method) 

- Pseudo-invariant targets 

- Ground-based atmosphere (aerosol) measurements 

- Scene-based aerosol estimates (based on dark pixels) 

- Climatological atmosphere 

- Satellite-based atmospheric measurements 



This presentation will focus on preliminary results of only the 
satellite-based atmospheric correction algorithm. All algorithms 
will be evaluated in the coming year. 
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Atmospheric Correction Prototype 

Algorithm 

Stennis Space Center 

• Leverage JACIE commercial imagery radiometric 
characterizations 







- IKONOS, QuickBird, OrbView-3 (future) 

Use daily coverage from MODIS to provide input data for 
atmospheric correction 


- MOD04 Aerosol Optical Thickness- 

- MOD05 Total Precipitable Water (Water Vapor) 

Generate MODIS-like products Depth 

- Surface Reflectance (MOD09) 



Gridded Vegetation Indices - NDVI (MODI 3) 


Water Vapor, 
Near Infrared 
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Atmospheric Correction Approach 

Stennis Space Center 
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MODIS 

(Moderate Resolution Imaging Spectroradiometer) 

^ ^ ^ ^ Stennis Space Center 



MODIS provides long-term observations from which an enhanced knowledge of 
global dynamics and processes occurring on the surface of the Earth and in the 
lower atmosphere can be derived. 


MISSIONS: 

• Terra - Dec 1999 

• Aqua - May 2002 

HERITAGE: 

• AVHRR 

• High Resolution Infrared 
Radiation Sounder (HIRS) 

• Landsat TM 

• Coastal Zone Color Scanner 

LINKS: 

• Sensor Site: 
http://modis.asfc.nasa.gov/ 


PRODUCT SUMMARY: 

• Congruent observations of high-priority v 
atmospheric, oceanic, and land-surface 
features 

VITAL FACTS: - ' 

• Instrument: Whiskbroom imaging radiometer 

• Bands: 36 from 0.4 and 14.5 pm 

• Spatial Resolution: 250 m (2), 500 m (5), 
1000 m (29) 

• Swath: 2,300 km (±55°) from 705 km 

• Repeat Time: Global coverage in 1 to 2 days 

• Design Life: 6 years 


• Data Sites: 

http://daac.qsfc.nasa.gov/ (ocean and atmospheric) 
http://edcdaac.usqs.gov/main.html (land) 


OWNER: 

• U.S., NASA 
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Spherical Albedo Formulation 



Stennis Space Center 


The spherical albedo approach approximates the signal observed by 
the satellite as the summation of the components illustrated below 



Solar Irradiance 

<( )> 


Atmosphere with 
spherical albedo, s 
(aerosols, 
molecules, 
pressure, 

temperature, 

humidity) Background 


Satellite 


B PbgSPb 9 


Pbg 

MgtSW 


Target 

L TO A ~ K + Ap m + Ap sp b + Ap s 2 p b 2 + ... + Bp bs + Bp b sp b + Bp,s 2 p 2 + 
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tmospheric Correction Approximations 

Stennis Space Center 


Spherical Albedo formulation simplifies to: 


Ltoa — A > ~*~ 


Mg, B P, 


+ 


bg 


1 - sp hg 1 -sp bi 



Knowledge of atmosphere 
and adjacency 

Knowledge of atmosphere 


Where : 

p = Target reflectance 

p bg = Background reflectance 

L toa =Top of atmosphere (at - sensor) radiance 

A, B, s, and L 0 are constants that depend on 

atmospheric properties and geometry 
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Adjacency Effects 

Stennis Space Center 

• Adjacency effects are caused by complicated multiple 
scattering in the atmosphere-land surface interactions 

- Dark pixels appear brighter and bright pixels appear darker 

- Significant in turbid atmospheres over highly heterogeneous 
landscapes 

• Different methods have been employed for removing this 
effect 

- Atmospheric point spread function-PSF (Environmental 
Function) 

- Empirical formula 
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Spherical Albedo Benefits 


Stennis Space Center 

• Commonly used and found throughout the literature 

• Allows for analytical determination of target 

albedo/reflectance values v 

• Computationally efficient 
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Atmospherically Corrected Imagery 



Stennis Space Center 


Stennis Space Center 
January 15, 2002 



Includes material © Space Imaging, LLC 


IKONOS CIR image 
(rgb=431) 
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Atmospherically Corrected NDVI from 
IKONOS Imagery 


Stennis Space Center 



Stennis Space Center 
January 15, 2002 


Includes material © Space Imaging, LLC 
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NDVI from IKONOS Imagery 



Includes material © Space Imaging, LLC 


Stennis Space Center 
January 15, 2002 


Stennis Space Center 
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Count 


NDVI Histogram Comparison 



Stennis Space Center 

1000 r 1 



NDVI 
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Stennis Space Center 


Atmospheric Correction Prototype 
Algorithm Verification 



Scene Selection for Atmospheric 
Correction Algorithm Verification 

SteHIliS SpdCe Center 


• Criteria 
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n 
h 
E 
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NIR Band Calibration Summary - Average Spectral Radiance 


SSC,Stennis,52%U»rp. 1 /1 0/04 
SSC,Stannls,34%tarp, 1 /1 0/04 
SSC, Stanmc.22%tarp. 1 /1 0/04 
SSC.Stennis.3.5%twp, 1/1 0/04 
SSC.Stenms.Concrote. 1/1 0/04 
SSC .Stonms, Grass, 1 /1 0/04 
SSC, Brookings, Grass, 8/2 3/03 
SSC . Brookings. 52%t<jr p, 9/1 5/03 
SSC. Brooklngs.3 5%tarp,9/1 5/03 
SSC,Brookingc,Grace.fl/1 5/03 
SSC, Brookings, Grass, 10/21/03 
UofA.WSMR Playa.8/22/03 
UofA, I vanpah, Playa, 1 2/1 5/03 
UofA,RRV,P1ayo, 1/7/04 
DG Cal Curve 
SSC Cal Cun/e 
SSC Cal Curve ± 1o 


Available ground-truth reflectance and 
measurements 

- Available radiometric calibration 
coefficients 

- MODIS overpass close in time to 
IKONOS/QuickBird overpass 


atmospheric 
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Selected Scene Matrix 


Stennis Space Center 


Location 

Date 

HSR 

Sensor 

Sensor 

Az/EI 

Ground Truth 

SSC, MS 

Jan. 15, 
2002 

IKONOS 

113.0/77.2 

deg 

Targets = 3 tarps (3.5, 22, 52), grass, concrete 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR, 
pressure, radiosonde, BRDF 

SSC, MS 

Feb. 17, 
2002 

IKONOS 

100.7/81.9 

deg 

Targets = 3 tarps (3.5, 22, 52), grass, concrete 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR, 
radiosonde, BRDF 

SSC, MS 

Feb. 17, 
2002 

QuickBird 

10.5/67.3 

deg 

Targets = 3 tarps (3.5, 22, 52), grass, concrete 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR, 
radiosonde, BRDF 

Brookings, SD 

July 20, 
2002 

QuickBird 

349.8/64.1 

deg 

Targets = 2 tarps (3.5, 52), grass 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR, 

radiosonde 

Brookings, SD 

Sept. 7, 
2002 

QuickBird 

191.0/74.9 

deg 

Targets = 2 tarps (3.5, 52), grass 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR 

SSC, MS 

Nov. 14, 
2002 

QuickBird 

274.8 / 79.4 
deg 

Targets = 3 tarps (3.5, 22, 52), grass, concrete 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR, 
pressure, radiosonde, BRDF 

Brookings, SD 

Aug. 23, 
2003 

QuickBird 

148.3/76.8 

deg 

Targets = grass 

ASD FieldSpec FR Spectroradiometer, ASR/MFRSR 

Brookings, SD 

Oct. 21, 
2003 

QuickBird 

279.5/81.3 

deg 

Targets = grass 

ASD FieldSpec FR Spectroradiometer, ASR/MFRSR 

SSC, MS 

1 L/. Kjyj'T 

Jan. 10, 
2004 

QuickBird 

230.7 / 89.2 
deg 

Targets = 4 tarps (3.5, 22, 34, 52), grass, concrete 
ASD FieldSpec FR Spectroradiometer, ASR/MFRSR, 
pressure, radiosonde 


Satellite Overpass Times 


Stennis Space Center 


Location 

Date 

HSR 

Sensor 

HSR Satellite 
Overpass Time 

MODIS 

Overpass Time 

SSC, MS 

Jan. 15,2002 

IKONOS 

16:44 

17:10 

SSC, MS 

Feb. 17, 2002 

IKONOS 

16:47 v 

16:14 

SSC, MS 

Feb. 17, 2002 

QuickBird 

16:45 

16:14 

Brookings, SD 

July 20, 2002 

QuickBird 

17:26 

17:42 

Brookings, SD 

Sept. 7, 2002 

QuickBird 

17:22 - 

16:47 

SSC, MS 

Nov. 14, 2002 

QuickBird 

16:44 

16:25 

Brookings, SD 

Aug. 23, 2003 

QuickBird 

17:07 

16:57 

Brookings, SD 

Oct. 21,2003 

QuickBird 

17:11 

18:17 

SSC, MS 

Jan. 10,2004 

QuickBird 

16:30 

17:27 


November 10, 2004 


21 



Algorithm Validation using NASA/JACIE 
Generated Ground Truthing Datasets 

Steniiis Space Center 


Laboratory Measurements 


ASD FieldSpec FR Spectroradiometer calibrations 
BRDF laboratory measurements 


• Field Measurements 

- Radiometric calibration tarps, grass, and concrete targets 

- In-field calibrated sun photometers 

In-field setup to check atmospheric model parameters 
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Calibration and Characterization of 
ASP FieldSpec Spectroradiome ters 


Stennis Space Center 


• NASA SSC maintains four ASD 
FieldSpec FR spectroradiometers 

Laboratory transfer radiometers 

Ground surface reflectance for V&V 
field collection activities 

• Radiometric Calibration 

NIST-calibrated integrating sphere- 
serves as source with known spectral 
radiance 

• Spectral Calibration 

Laser and pen lamp illumination of 
integrating sphere 

• Environmental Testing 

Temperature stability tests performed 
in environmental chamber 
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Laboratory BRDF Measurements 


Stennis Space Center 


• Purpose 

- Laboratory BRDF measurements 
are used to correct ground-based 
reflectance measurements for 
satellite viewing and for solar 
illumination geometry 

• Method 

- Collimated FEL lamp source 

- NIST-calibrated Spectralon® 
panel serves as reference 

- Goniometer-mounted sample 
controls illumination geometry 

- Optronics OL750 hyperspectral 
instrument measures spectra 
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Algorithm Verification Case Matrix 

Stennis Space Center 

• Three different atmospheric correction approximations 

- Case (1 ) Planetary reflectance 

- Case (2) Spherical albedo w/knowledge of atmosphere 

- Case (3) Spherical albedo w/knowledge of atmosphere & 

adjacency 

• Three different sets of data used as input into 
approximation 

- Case (a) ground based-sun photometer (aerosol), TOMS 

(ozone), Radiosonde (water vapor)- 

- Case (b) MOD04 (aerosol), TOMS (ozone), Radiosonde (water 

vapor) 

- Case (c) MOD04 (aerosol), MOD05 (water vapor), TOMS 

(ozone) Operational Case 

• Nine different scenes 

9 cases (1) + 17 cases (2b/2c) + 20 cases (3a/3b/3c) - 46 cases 
Planetary Spherical Albedo Spherical Albedo 

November 10, 2004 W/adjdCeilCy 25 



Optical Depth 



MODIS/Sun Photometer Comparison 

^ Stennis Space Center 


Total Optical Depth 
Stennis Space Center 
January 10, 2004 
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Reflectance Map Over SSC Tarps 

^ Stennis Space Center 


January 10, 2004 - Case 3c ( Operational input of best approximation) 


Blue Band Reflectance 


0.2 


0.4 
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Z 3 


0.6 


0.8 


0.2 0.4 0.6 0.8 
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Green Band Reflectance 




NIR Band Reflectance 













Measured and Calculated Reflectance 
Values of 52% Tarp 


Algorithm 

Approximation Effects 


Stennis Space Center 

Stennis Space Center 
January 10, 2004 


CASE 

reflectance 

(Case reflectance - meas reflectance) 

BLUE 

GREEN 

RED 

NIR 

RMS 

Measured Reflectance 
(Truth) 

0.52 

0.52 

% 

0.53 

0.53 

— 

Case 1 

(Planetary reflectance) 

0.46 

(-0.06) 

0.46 

(-0.06) 

0.48 

(-0.05) 

0.50 

(-0.03) 

0.05 

Case 2c 

(Operational - no adjacency) 

0.41 

(-0.11) 

0.44 

(-0.08) 

0.47 

(-0.06) 

0.50 

(-0.03) 

0.07 

Case 3c 

(Operational - w/adjacency) 

0.50 

(-0.02) 

0.50 

(-0.02) 

0.51 

(-0.02) 

0.52 

(-0.01) 

0.02 


Important to take into account the adjacency effect 
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Measured and Calculated Reflectance 

Values of Grass Target 

— Stennis Space Center 


CASE 

reflectance 

(Case reflectance - meas reflectance) 

BLUE 

GREEN 

RED 

NIR 

RMS 

Measured Reflectance 
(Truth) 

0.05 

0.07 

% 

0.10 

0.18 

— 

Case 1 

(Planetary reflectance) 

0.15 

(0.10) 

0.13 

(0.06) 

0.13 

(0.03) 

0.20 

(0.02) 

0.06 

Case 2c 

(Operational - no adjacency) 

0.06 

(0.01) 

0.08 

(0.01) 

0.11 

(0.01) 

0.20 

(0.02) 

0.01 

Case 3c 

(Operational - w/adjacency) 

0.06 

(0.01) 

0.08 

(0.01) 

0.11 

(0.01) 

0.20 

(0.02) 

0.01 


Stennis Space Center 
January 10, 2004 


Algorithm 

Approximation Effects 
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Measured and Calculated Reflectance 

Values of 52% Tarp 

— — Stennis Space Center 


MODIS Input 
Parameter Effects 


Stennis Space Center 
January 10, 2004 


CASE 

reflectance 

(Case reflectance - meas reflectance) 

BLUE 

GREEN 

RED 

NIR 

RMS 

Measured Reflectance 
(Truth) 

0.52 

0.52 

% 

0.53 

0.53 

— 

Case 3a 

(Meas. aerosol and water) 

0.52 

(0.00) 

0.52 

(0.00) 

0.53 

(0.00) 

0.53 

(0.00) 

0.00 

Case 3b 

(MOD04 & meas. water) 

0.50 

(-0.02) 

0.50 

(-0.02) 

0.51 

(-0.02) 

0.52 

(-0.01) 

0.02 

Case 3c 

(MOD04 & MOD05) 

0.50 

(-0.02) 

0.50 

(-0.02) 

0.51 

(-0.02) 

0.52 

(-0.01) 

0.02 
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Measured and Calculated Reflectance 
Values using Operational Inputs 

Stermis Space Center 


MOD04 - aerosol 
MOD05 - water vapor 


Stennis Space Center 
January 10, 2004 


CASE 3c 

reflectance 

(Case reflectance - meas reflectance) 

BLUE 

GREEN 

RED 

NIR 

RMS 

3.5% Tarp 

(0.04, 0.04, 0.03, 0.03) 

0.05 

(0.01) 

0.06 

.(0.02) 

0.05 

(0.02) 

0.06 

(0.03) 

0.02 

22% Tarp 

(0.22, 0.22, 0.21, 0.20) 

0.24 

(0.02) 

0.24 

(0.02) 

0.23 

(0.02) 

0.23 

(0.03) 

0.02 

34% Tarp 

(0.31, 0.31, 0.31, 0.30) 

0.34 

(0.03) 

0.34 

(0.03) 

0.34 

(0.03) 

0.34 

(0.04) 

0.03 

52% Tarp 

(0.52, 0.52, 0.53, 0.53) 

0.50 

(-0.02) 

0.50 

(-0.02) 

0.51 

(-0.02) 

0.52 

(-0.01) 

0.02 

Grass 

(0.05, 0.07, 0.10, 0.18) 

0.06 

(0.01) 

0.08 

(0.01) 

0.11 

(0.01) 

0.20 

(0.02) 

0.01 

Concrete 

(0.11, 0.13, 0.16, 0.19) 

0.12 

(0.01) 

0.15 

(0.02) 

0.17 

(0.01) 

0.21 

(0.02) 

0.01 
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Operational Algorithm Verification 

Summary (Case 3c) 

Stennis Space Center 


Location 

Date 

Sensor 

Average RMS for 
all targets 

SSC, MS 

Jan. 15, 2002 

IKONOS 

0.02 

SSC, MS 

Feb. 17, 2002 

IKONOS v 

0.04 

SSC, MS 

Feb. 17, 2002 

QuickBird 

• /* % 

0.01 

Brookings, SD 

July 20, 2002 

QuickBird 

0.01 

Brookings, SD 

Sept. 7, 2002 

QuickBird 

0.02 

Brookings, SD 

Aug. 23, 2003 

QuickBird 

0.00 

Brookings, SD 

Oct. 21, 2003 

QuickBird 

0.02 

SSC, MS 

Jan. 10, 2004 

QuickBird 

0.02 
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Summary/Future Direction 



Stennis Space Center 


• MODIS products (MOD04, MOD05) provide the 
necessary inputs to generate high-spatial-resolution 
reflectance products under many conditions 

- Average RMS differences range between 0.00-0.04 for the eight 
datasets evaluated (Case 3c-Operational input to best 
approximation) 

« v * 

• Adjacency can be an important component that needs to 
be accounted for to minimize errors 

• Future Activities/Recommendations 

- Evaluate alternate algorithms 

- Compare algorithm results to MODIS products (MOD09, 

MOD13) 

- Compare algorithm results to commercial atmospheric correction 
algorithm results (FLAASH, ACORN, ATCOR ...) 
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